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I. INTRODUCTION 



Recently the problem of dynamics of nonlinear lattices under strong and rapid modula- 
tions of parameters has attracted a lot of attention. Two systems have been analyzed. The 
first one is the diffraction-managed array of optical waveguides, with diffraction varying pe- 
riodically along the beam propagation pj. The model is described by the discrete nonlinear 
Schrodinger (DNSL) equation with rapidly and strongly varying in time tunnel-coupling 
between sites coefficient c(t): 

iUnt + ^c(^)(u n+ i + U n -i) + 2\u n \ 2 u n = 0, (1) 

where e<l. The analysis exhibits the existence of a new type of discrete spatial optical 
solitons with beam width and peak amplitude evolving periodically during propagation. 

The second system is the Bose-Einstein condensate in a periodic (in space) potential 
with a varying (in time) scattering length. In the tight-binding approximation this system 
is described by the DNLS equation [2] with strongly and rapidly varying in time nonlinearity 
coefficient n(t): 

iu nt + {u n+ i + u n _i) + -n(-)\u n \ 2 u n = 0. (2) 

It ™ show, that this S yste m supports nonBn^Uanased d.screte soBtons 3 . In a more 
general context it is of interest to investigate the influence of rapid perturbations on the 
dynamics of discrete solitons in nonlinear lattices. The case of strongly and rapidly varying 
external drivers is particularly important for applications. This problem is encountered 
both in the study of the dynamics of a magnetic flux quantum in an array of long Josephson 
junctions with varying ac current [4| and in the evolution of an optical field in a nonlinear 
chain of resonators or microcavities in presence of pumping 5, 6, 7, sj]. 

In this paper we consider the influence of a rapid strong drive on discrete bright solitons 
and cnoidal waves of the Ablowitz-Ladik (AL) and the DNLS equations with damping. 
Although the AL system has scarce physical applications it has many advantages from the 
analytical point view such as the complete integrability of the unperturbed system, existence 
of moving discrete solitons etc. In some regions of the parameter space, the DNLS equation 
can be described as a perturbation of the AL model, a feature that we shall take advantage 
of in the following. The dynamics of discrete solitons in AL and DNLS equations under 

nn 

the influence of damping and slowly varying driving field has been studied in Refs. |9|, H0[. 
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Recently the influence of parametric drivers on the stability of strongly localized modes 



of the DNLS equation near the anti-continuum limit has been investigated in 



The 



stability of solitons in the continuous parametrically driven NLS equation has been studied 



m 
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131 ] . Here we address a general discrete nonlinear system with a strong rapid drive 
modeled by the following equation: 

iu nt + [U n+ i - (2 +u)u n ] + |«n| 2 [(l - x) + ^n-l) + 2 X«n] = ~/(^) ~ '7%) ( 3 ) 

where / is a zero-mean periodic function with period 1 that describes the rapid drive, and 
the small parameter e is the period of the drive. Here the parameter 7 > denotes the 
damping term, uo is the propagation constant in optics (chemical potential in BECs), and 
the parameter x £ [0, 1] characterizes the type of nonlinearity. For \ = the nonlinearity 
is of intersite type as in the AL model while for x = 1 we have the onsite nonlinearity of 
the DNLS equation. In absence of strong rapid drive and damping Eq. (j3J) coincides with 
the Salerno model [yj] which interpolates between the AL model and the DNLS model. We 
show that by averaging out the fast time scale one can reduce Eq. ([3]), for particular choices 
of parameters, to the parametrically driven AL and DNLS equations with damping. The 
existence of parametric discrete bright solitons and cnoidal waves of these equations is then 
investigated and the stability properties is analyzed both analytically and numerically. We 
find that the analytical predictions obtained from the averaged equations by means of a 
perturbation scheme based on the inverse scattering method are in good agreement with 
direct numerical simulations of the problem with rapidly varying drive and damping. 

Finally, we remark that the physical systems where parametric discrete solitons of the type 
considered in this paper can be realized are chains of linearly coupled nonlinear microcavities 



|15j and nonlinear waveguide arrays with dielectric mirror at the ends, driven by an external 
time-dependent field {7, 8]. In this case the equation describing the discrete cavity solitons 
has the form of the parametric DNLSE 

iu nt + u)u n + ryw n + a\u n \ 2 u n + C(u n+1 + w n _i - 2u n ) = F n (t), (4) 

where u is the detuning from the linear resonance parameter, 7 is the effective damping 
parameter, C is the effective coupling between adjacent waveguides, F n (t) is the input field 
in the n-th waveguide. Injecting a homogeneous in space and rapidly varying in time field, 
we can generate discrete parametric soliton in this system. 



The paper is organized as follows. In Section II we derive the averaged DNLS-type 
equation for the strongly and rapidly varying external drive model. In Section III we analyze 
the solitons in a damped AL system with parametric drive. We use the perturbation theory 
based on the Inverse Scattering Transform (1ST) and study the stability region of parametric 
discrete solitons. Periodic solutions of the damped AL equation with parametric drive are 
also found in this section. The discrete soliton dynamics in parametrically driven DNLSE 
is investigated in Section IV. 



II. AVERAGING 

We look for the solution of Eq. ([3]) in the form 

Un (t) = u M(t, t -)+euM(t,~) + ... (5) 

where uh ^ i Un \ ••• are periodic in the argument r = t/e. We substitute this ansatz into 
Eq. ([3]) and collect the terms with the same powers of e. We obtain the hierarchy of equations: 

iuf T = m 

+ [«Sa + «£i - (2 + ^R 0) ] + k 0) H(l - + 41) + 2rf] = -iiu^ - iu« T 

The first equation imposes the form of the leading-order term 

u^(t,r) = -iF{r) + a n (t) , 

where F(r) = f(s)ds is the antiderivative of / and a n depends only on the slow variable 
t. The second equation is the compatibility equation for the existence of the expansion ([51). 
By integrating over a period in r, we obtain: 



ia nt + [a n+1 + an-! - (2 + u)a n ] + (1 - x) (K - iF(-)\ 2 (a n+1 + a n _i - 2iF( 

+2 X (\a n - iF{-)\ 2 {a n - iF{-))) = -i 7 a n 

which gives 

ia n t + [a n +i + a n _i - (2 + u)a n ] + K| 2 [(l - x)( a n+i + a n -i) + ^X a n] 

= -ija n + 8an~ - -[(1 - x)( a n+i + a-n-i) + 2(1 + X ) a n), (6) 

where 5 = 2(F 2 ). 



Thus, the averaging method applied to Eq. ([3]) leads to a parametrically driven nonlinear 
lattice equation which reduces to the Salerno model jjjj in absence of perturbations. It is 
possible to consider a random forcing instead than a periodic one. More precisely, we get the 
same result if the source /(r) is a colored noise with coherence time of order one. However, 
the power spectral density of the source should vanish at zero-frequency. Otherwise it would 
appear a phase diffusion and this would destroy the stability of the stationary solution that 
we will introduce next. 

We introduce the rescaled time T — t(l + 5(1 — x)/2) and the rescaled function A n = 



a n j a/1 + <5(1 — x)/2. The averaged equation for the function A n (T) has the form 

%A nT + [A n+1 + A„_! - (2 + n)A n ] + \A n \ 2 (A n+1 + A n ^) = R n , (7) 

where 

a= J , r= „ 7 , , n u ~ 26 



l + «J(l-x)/2' l + 5(l- X )/2' l + 5(l- X )/2' 

and the term in the right-hand side is given by 

R n = -iVA n + AA^ + x \A n \ 2 (A n+1 + A n ^ - 2A n ). (8) 

The left hand side of ((7j) is the Ablowitz-Ladik (AL) equation, that is completely integrable. 
The right-hand side can be seen as a perturbation of this system. We shall use the perturbed 
Inverse Scattering Transform (1ST) to study the evolution dynamics of AL solitons driven 
by the perturbation R n . 

III. THE DAMPED AL SYSTEM WITH PARAMETRIC DRIVE 

We consider in this section the case \ = 0, that is the AL model with damping and rapid 
drive. Therefore, we consider the perturbed AL equation (j7]) with the perturbation 

Rn = -iTA n + AA^. (9) 
A. Perturbed Inverse Scattering Transform 

We assume that the damping parameter 7 and the parametric drive parameter S are small 
(but uj can be of order one). Therefore, T and A are small, and, following 

0,0,0, the 



5 



evolution equations for the soliton parameters in the adiabatic approximation have the form 

sinh/3 . sinh/9 *^ (n — x) coshpYn — x)lm(r n ) 

x T = 2 — - — sin cH — > 

R R t-^ 



P P n^oo cosh/5(n + 1 — x) cosh(/9(n — 1—x) 

(3 - sinh/3 cosh/j(n - x)Im(r ra ) 

^ cosh/3(n + 1 — x) cosh/3(n — 1 — x) ' 



n=— oo 

oo 



. , a v-^ smh/3(n - x)Re(r n ) 

a T = smhp > — — - r — - , 

cosh pin + 1 — x) coshpln — 1 — x) 

n=— oo v ' 

sinh/3 (n — x) sinh/3(n — x)Re(r n ) 

o T = 2 cosh p cos a + 2 — - — asm a — 2 — J 2 + smhp > — — — — — 

p ^-^ cosh p(n + 1 — x coshpln — 1 — x) 

n=—oo v / < \ / 

-cosh/? V cosh/3(n-a:)Re(r TO ) 

^ cosh /3(n + 1 — x) cosh /3(n — 1 — x) 



n=— oo 
oo 



sinh/3 ^-^ (n — x) cosh/3(n — x)Im(r n ) 



/3 f-^ cosh (3(n + 1 — x) cosh /3(n — 1 — x) 



where r n = R n exp(—ia(n — x) — id). Using standard analytical tools (Poisson summation 
formula and residue theorem), we can compute the right-hand sides of these equations. The 
equation for (5 takes the form: 

fa = APW + rG (/3) (10) 



where 



with 



P {/3) = -sinh 2 /5 //3(2a + 27rs)sin(27rsx + 2cr) , 

s=— oo 

G (/3) = -2tanh/3 



27rsina 



/3sinh(2/3) sinh(|f)' 
The equation for a has the form: 

a T = AP (a) + TG {a) (11) 

where 

oo 

P (a) = - sinh 2 /3 ^ ify^a + 2vrs) sin(27rsx + 2a) , 

S = — OD 

G {a) = , 



6 



with 

2tt sin 2 (a/2) 



K p {a) 



/3sinh'/3sinh(§) 
The equation for the soliton center x has the form: 

XT = 2 ^l sina + Ap{x) + TG{x) 



where 



p{x) = _swhj_ j /3 (2a + 2ns)cos{2nsx + 2a) , 



s=—oo 

2 



with Jp (a) = —I 'J a) and 



(n — x) An ^-^ sin(27rsa;) 

~ n ^ cosh(/3(n + 1 - rr) cosh(/3(n - 1 - x) ~ /3sinh(2/3) ^ sinh(^) 

Finally, the equation for the soliton phase a has the form: 

a T = 2 cosh /3 cos a + 2^^-a sin a - 2 - ft + AP (,t) + rG ((T) 



where 



P (ff) = -sinh/3 L /3 (2a + 2tts) cos(a + ws + 2a + 2nsx) 

S = — QO 
OO 

+ sinh/3 ^ K' fS (2a + 2irs)sm(2irsx + 2a) 

S— — OO 

OO 

-sinh/? ^(2a + 27rs) cos(27rsx + 2(j) 

s=— OO 

a sinh 2 /3 *^ w . n . . _ „ . 

2^ Jj(2a + 2tts) sin(27rsa: + 2a) , 



2/5 



s=— OO 

2 



sinh 8 , . . 
G = -a—^-^x) , 



with 



. . 27rsin('a/2 v ) 



B. Parametrically driven AL solitons 



The system ( ITOlfTB"]) has two fixed points if T < A (which is equivalent to 7 < 5) and 
Q + A > (which is equivalent to uj — 5 > 0): 



a± = 0, sin(2a ± ) = -^-, cos(2a ± ) = 1 - ^ , (14) 
(3± = arccosh ^1 + — + — cos(2cr ± )^ . (15) 

In fact, the second point (with subscript — ) exists only if Q + Acos(2cr_) > 0. The center 
of the soliton can be arbitrary. Note that these values correspond to a fixed point of the 
equation (J7|). The soliton ± is of the form 

AJT) = sinh/3 ± c , CT+ 
n cosh [f3±(n — x±)} 

We can investigate the linear stability of these solutions. The linearization of the nonlinear 
system ( TT0lfT3T) around the parameters of the stationary solitons gives the linear system 



Pi T = — 2 A sinh 2 P±ip/3 ± (x±) cos(2cx±)a;i — 4A tanh/3± cos(2cr±)cri , (16) 
a>i T = — 2Tai , (17) 

a 1T = 2sinh(3 ± (3 l -2Ta 1 , (18) 

sinh/3± A sinh 2 /3 ± tt 2 + A(3 2 ± 

xit = 2—5 ai-A — 2 . . cos(2o- ± )a! 

p± p± 3p± smh(2p ± ) 

sinh 2 /3 ± J3± sinh(2/? ± ) - sinh 2 /3 ± 

-r — - — i) P± (x±)xi - r — 2 ip{3 ± {x ± )(3 1 . (19) 

p± p± 

The equation for x\ is decoupled from the other ones. The three eigenvalues for the 3x3 
system for (ax, Pi, ax) are 

Ai 1} = -2T, 

A^ 2) = -r + - 8A sinh 0± tanh (3± cos(2(T ± ) , 
\® = -r - ^r 2 - 8A sinh (3± tanh (3± cos(2(j ± ) . 

The fixed point is stable if the real parts of the eigenvalues are nonpositive. This shows that 

(2) 

the fixed point labeled — is not stable, since A_ > 0, while the fixed point labeled + is 
stable. 
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FIG. 1: Here x = 0, u = 1, 5 = 0.022, and 7 = (i.e. A = 0.0218 and T = 0). The initial 
condition is a soliton with a = a + , x + = 0, a = 0, and /? = /?+ + 0.02 ( (3 + = 0.948). We plot the 
solution phase obtained from the numerical integration of Eq. (|6|) (dashed line) and compare with 
the theoretical oscillation obtained from (|16M18p (solid line). The theoretical oscillation period (in 
the t-variable) is 2vr/ft + /(l + 5/2) = 16.6. 

In absence of damping T = (which is equivalent to 7 = 0), the soliton parameters (3 
and a (amplitude and phase) oscillate with the frequency (in the T-variable) 

Q + = y/8 A sinh ( 5 + tanh/3+ (20) 

while the parameter a (velocity) is constant. This means that the stationary soliton is 
stable. This also shows that stable slowly moving breathers can propagate in the presence 
of parametric drive. 

We have performed numerical simulations of Eq. ([6]) to confirm these predictions. In 
Fig. [1] we consider a perturbation of the initial amplitude. The periodic oscillations of the 
soliton parameters have the predicted period (!20j) . In Fig. [2] we consider a soliton with a 
positive velocity. As predicted by the theory, this soliton can propagate in a stable way. 

In presence of damping the soliton parameters (3 and a oscillate with the frequency 

Q+ = ^/Wa 2 - r 2 sinh/3 + tanh/? + - T 2 . (21) 

These oscillations decay exponentially with the rate T. Besides, comparing (1201 and (12 lj) 
shows that the damping enhances the oscillation period. The soliton parameter a decays 
exponentially at the rate 2T. Therefore, the stationary soliton is very stable. However, the 
propagation of moving solitons is not supported, as the soliton velocity decays exponentially 
to 0. If we denote by a the initial value of the parameter a, then the input soliton converges 
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FIG. 2: Here x = 0, u; = 1, <5 = 0.022, 7 = 0. The initial condition is a soliton with a = <r+, 
x+ = 0, a = 0.1, and f3 = /3+. We plot the soliton profiles |a n (£)| at different times, which exhibits 
the stable propagation of the moving soliton. 




200 



FIG. 3: Here x = 0, u = 1, 5 = 0.022, 7 = 0.02. The initial condition is a soliton with a = <7+, 
x + = 0, a = 0, and /3 = f3 + + 0.02 (/?+ = 0.942). We plot the solution phase obtained from 
the numerical integration of Eq. ([6]) (dashed line). The observed oscillations and damping are 
correctly predicted by the model (|16H18p (solid line). The period is 27r/S7 + / (1 + 5/2) = 25.9 and 
the exponential decay rate is 7 = 0.02 (solid line). 



to the stationary form 



with 



A r . 



sinh /3 4 



cosh [P+(n — x + — xp)\ 
sinh (3 + sin ao 



3 t(T+ 



Xp 



(22) 



In Fig. [31 we plot the damping of the oscillations of the soliton parameters and compare 
it to the theoretical formula. In Fig. H] the trapping of a soliton with an initial velocity is 
shown. The final position of the soliton is given by f[2"2"j) . 

We have also simulated the original equation ([7j) with the external drive /(r) = sin(27rr) 
and e = 0.125. We plot one of the obtained results in Fig. [5] (to be compared with Fig. [3]), 
which shows full agreement. 
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FIG. 4: Here % = 0, u = 1, 8 = 0.022, 7 = 0.02. The initial condition is a soliton with a = a+, 
x+ = 0, a = 0.1, and /? = /?+. We plot the soliton profiles |a n (i)| at different times, which 
exhibits the trapping of the moving soliton. The solid line is the theoretical stable stationary 
soliton centered at xf = 5.76 given by (f22j) . 
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FIG. 5: Numerical simulation of Eq. ©. Here x = 0,u = l,5 = 0.022, 7 = 0.02, /(r) = sin(27rr) 
and e = 0.125. The initial condition is a soliton with a = <r+, x + = 0, a = 0, and /? = /?+ + 0.02. 
We plot the phase of ao(t) := uo(t) + isin(167rt) (dashed line). The observed oscillations and 
damping are correctly predicted by the model, and correspond exactly to the ones of the solution 
of the averaged equation (see Fig. [3]). 

C. Periodic solutions of the damped AL equation with parametric drive 

In the following we discuss exact periodic solutions of the parametrically damped and 
driven AL equation 

iA nT +[A n+1 -(2 + Q)A n + A n _ 1 ] + [i v +{l+i v )\A n \ 2 }{A n+1 + A n _ 1 ) = -iTA n + AA^,, (23) 

where rj > models a small dispersive and nonlinear damping. To this regard we consider 
periodic stationary solutions of the form 

A n = Ae ia dn [f3(n + x F ), m] . (24) 
11 
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FIG. 6: Left panel: Square modulus of the dn solution in Eq. (I24p with m = 0.1, j3 = 0.161244 
(Np = 20), xp = 0, on a line of 160 points. The parameters are fixed as A = 0.04, V = 0.02, 
£1 = 0.01558, rj = 0. Right panel: Square modulus of the cn solution in Eq. (|28|) with parameters 
fixed as in the left panel. 

Direct substitution of Eq. fl24l into Eq. fl23l) shows that a solution is indeed obtained 
provided the two following relations are satisfied by the soliton parameter f3 and the modulus 
m: 

2" 



„ dn( B,m) 
2 + VL- 2 — — — - + A 

cn(/3, m) 2 

(3N P = 2K(m), 



{ _ , £ 2 r] dn(/3, m) 



A Acn(/3,m) 2 



1 1/2 



0. 



(25) 
(26) 



where N p is the number of sites in one period and K(m) is the complete elliptic integral 
of the first kind. Here m G (0, 1) and iVp must be a positive integer, so that there exists 
a numerable set of pairs (/3,m) that satisfy the conditions ()25H26p . If the conditions are 
fulfilled, then the periodic function (I2"4~j) is a solution of (I2"5j) with the amplitude A and 
phase a given by 

sn(/5, m) 



A 



1 (Y 2r?dn(/3,m) 
o — — arcsm — + — -- 

2 \A Acn(/3,m) 2 



(27) 



cn(/3, m) 

while the soliton center xf is arbitrary. 

Another periodic solution of Eq. (|23l) can be constructed by replacing the function dn 
into the ansatz (T2"4~|) with the elliptic cosine cn: 



A n = Ae ia cn [f3{n + x F ),m] . 
In this case the two conditions to be satisfied by (3 and m are 



2 + ^-2 



cn(/3, m) 
dn(/3, m) 2 

/5A^ P = 4FT(m), 



A 



r 2rjcn(f3 , m) 



A Adn(/3,m) s 



1/2 



(28) 

(29) 
(30) 
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and the soliton amplitude and phase are 

A = ±^n—— -, a- 

dn(p, m) 



1 / T 2/7 cn(/3, m) 
•- arcsm — + -— - — — r— 

2 \A Adn(/3,m) 2 



(31) 




20 40 60 80 100 120 140 160 




200 E 



40 80 120 160 

n 



FIG. 7: Time evolution of \A n \ 2 for of the dn (left panel) and cn (right panel) solutions in Fig. [6]. 



It may be worth noting that in deriving these solutions usage of the following identities 
for the Jacobi elliptic functions [3] have been made: 

2dn(a, m) dn(x, m) 



dn(a; + a, m) + dn(:r — a, m) 
cn(x + a, m) + cn(x — a, m) 



cn(a, m) 2 + sn(a, m) 2 dn(x, m) 2 ' 
2 cn(a, m) cn(x, m) 



(32) 
(33) 



dn(a, m) 2 + m sn(a, m) 2 cn(x, m) 2 

Also notice that in the limit of infinite period (i.e. m — > 1) the above solutions both reduce 
to the AL soliton 



A r 



with the soliton parameter satisfying 



sinh(/3) 



cosh [/3(n + ip)] 



2 + Q + [A 2 - (r + 277 cosh(/3)) 2 ] 1/2 - 2 cosh(/3) = 0. 



(34) 



Similar periodic solutions exist also for the unperturbed AL equation [20| and for generalized 
AL equations with arbitrarily high-order nonlinearities 21 1. 

In Fig. [6] we depict the waveforms of the above solutions on a line of 160 points for the 
case A = 0.04, T = 0.02, Q = 0.01558, rj = 0, and the solution parameters are m = 0.1, 
(3 = 2K(m)/N p = 0.161244, A?p = 20. We find, by direct numerical integrations of Eq. (j23j). 
that for small values of the damping and parametric driver amplitude these solutions remain 
stable under very long time evolution. This is shown in Fig. [7] where the time evolution 
obtained from direct numerical simulations of Eq. ff23|) is depicted for parameters A = 0.15, 
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FIG. 8: Time evolution of the unstable dn solution in Eq. (|24p obtained for A = 0.15, T = 0.02, 
Q = —0.098438, r] = 0, and the solution parameters are as in Fig. [6j m = 0.1, f3 = 2K(m)/N p = 
0.161244, N p = 20. Notice the change in the \A n \ 2 scale from one panel to another. 



T = 0.02, Q = —0.098438, rj = 0, which supports the same dn and cn solutions as in Fig. El 
By keeping fixed the damping constant and increasing the amplitude of the parametric 
driver in a certain range, we find that these solutions remain stable for very long time, while 
outside of this range instabilities quickly develop. The development of the instability for out 
of range parameter is investigated in Figs. [HE 
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FIG. 9: Time evolution of the unstable cn solution in Eq. f)28|) when the parameters are fixed as 
in Fig. [6l Notice the change in the \A n \ 2 scale passing from one panel to another. 



Notice that while for the dn solution the instability suddenly sets in without any apparent 
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pattern, the instability of the cn solution seems to follow a precise pattern. In particular, 
from Fig. [9] we see that before the instability fully develops at time t « 200 the cn solution 
bifurcates into a period three solution at t ~ 80 which remains stable for a long time. The 
presence of a small dispersive nonlinear damping (controlled by the parameter if) effectively 
increases the stability of both the period one and the period three solutions, as one can see 
from Fig. [TUJ The scenario behind the development of the instability of the cn solution is 
quite interesting and deserves more investigations. 
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FIG. 10: Time evolution of the unstable dn solution in Eq. (|24f) obtained for rj = 0.02. Other 
parameters are fixed as in Fig. [9j 

IV. THE DAMPED DNLS SYSTEM WITH PARAMETRIC DRIVE 

We consider in this section that the nonlinear term is not of the AL type but it is a 
mixture of the AL cubic intersite nonlinearity and the onsite DNLS cubic term. We shall 
analyze this system by considering such a nonlinearity as a perturbation of the cubic AL 
nonlinearity. 

A. Perturbation theory for the cubic nonlinearity 

In this subsection we consider the perturbed AL equation (j7j) when the perturbation is 
given by 

Rn = X \A n \ 2 (A n+1 + - 2A n ) (35) 

If x — 0, then Eq. (j7]) with the perturbation (135]) is the integrable AL equation. If \ = 1, 
then it is the standard DNLS equation. In the adiabatic approximation, this perturbation 
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has no effect on the first-order differential equation for the amplitude soliton parameter 8, 
but the evolution equations for the parameters x, a, and a have corrective terms 

sinh/? . sinh 2 /? . 

x T = 2(1- X ) — sm a + 2 X — —sin a (36) 

«t = X-fMz) (37) 

ot = 2 cosh/? cos a + 2%(1 — cos a) sinh/3tanh/? 

sinh 3 sinh 2 

+2(1 - xh sm a— ^ + 2 X a sin a ^ -2-Q+ 2 X Qp(x) (38) 



where 



^ / \ 87r 3 s 2 sinh 2 8 . . s . , 



sinh 2/3 sinh 2 /? 
2vr 2 s/? 2 sinh /? cosh /3 + [vr 4 s 2 cotanh(^) - 2vr 2 s/?] sinh 2 



Qp{x) = — 1 H — sinh/?tanh/? 



+2 V g 2 cos(27rsx)(40) 

^ /3 4 sinh(^) 

We should keep in mind that the adiabatic approximation is valid when the perturbation fl35|) 
is small, which is true if x is small and 8 is arbitrary, or if x is arbitrary and 8 is small. In 
the following, we shall only keep the term s = 1 in the sums ( J39H401) . to simplify the algebra, 
although the analysis could be carried out with the full expressions. This simplification is 
consistent with the adiabatic approximation. 



B. Parametrically driven DNLS solitons 

We now consider that the perturbation R n is given by (jSJ), that is the sum of the cubic 
perturbation (|35|) and the parametric drive with damping (Q. In these conditions, there are 
two fixed points if T < A (which is equivalent to 7 < 5) and Q + A > (which is equivalent 
to u) - 5 > 0): 



o± = 0, sim> ± ) = ~ , cos(2a ± ) = ±yi-^, (41) 
cosh(/3 ± ) - 1 + X Q fe (0) = 5 + ^v/A^P . (42) 

The center of the soliton x± must be an integer as soon as X > 0- This is a manifestation 



of the Peierls-Nabarro barrier 116 



22j. Note that the function 8 i-» cosh(/?) - 1 + xQp(9) 
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is a one-to-one increasing function from (0, oo) to (0, oo) for any x > 0. Therefore, the 
parameter f3± is uniquely determined. The picture is the same as in the perturbed AL case. 
The only difference is a renormalization of the amplitude parameter (3± given by (1421) instead 
of $15]). 

We next perform the linear stability analysis of the fixed points. The linearization of the 
system of ordinary differential equations for the soliton parameters around the stationary 
points gives: 

Pi T = -4Atanh/5 ± cos(2<j ± )ai , (43) 

aiT = - 2Tai + 1 ^ x Pll H £) Xl > (44) 

a lT = 2 [sinh/3 + xdpQp(O)} p=0± Pi - 2Ta, , (45) 

sinh/?± sinh 2 /3 ± A tanh/? ±7 r 2 + 4^ 

x 1T = 2(1 -x) — - 5 ai + 2x-T2 — -a x - A — -3 cos(2cx±)ai 

p± p± cosh p± p± 6 



tanh/?± 4vr 2 

-r— 52 . rr^T ^- ( 46 ) 

P± smh(^) 

This 4x4 linear system can be decomposed into two 2x2 linear systems, for (fii, a\) and 
for (ai,Xi), respectively. It is then easy to compute the eigenvalues. They are: 



a£> = 
aL 2) = 


-r + n± 
-r - Ot 


Ai 3) = 




1 + 


Ai 4) = 


-r| 





tanh /3 ± 2tt 2 

1 smh(£ 
tanh /3 ± 2tt 2 



± 



^ sinh(g) 
where the complex numbers Q± and (l± are given by 

ft| = r 2 - 8A tanh f3± cos(2a±) [sinh (3 + xdpQ^)]p=p ± , (47) 

n| = r 2 fi- ta t^ 27r2 2 ^ 

I sinh(f 2 




4 smli p ± ( smh/3± smh fi ± 

/^smh(^) V P± /%cosh/? ± 

. tanh/3 ± 7r 2 + 4/% \ 
-A— gf± g-^ cos(2a ± ) J (48) 
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FIG. 11: Here % = 1 , w = 0.3, 5 = 0.022, 7 = 0.02. The T and t scales coincide. The initial 
condition is a soliton with a = <t+, x+ = 0, a = 0, and j3 = /3+ — 0.02 = 0.43. We plot in the left 
picture the argument of ao{t) (dashed line). The observed oscillations and damping are correctly 
predicted by the model. The period is 2ir/0, + = 43.9 and the exponential decay rate is T = 0.02 
(solid line). Note also, in the right picture, that the stationary profile is not exactly the sech 
predicted by the AL theory, but a slightly deformed version. 

The eigenvalues A± and Aj? describe the growth rates of the perturbations of the ampli- 
tude parameter (3 and phase parameter a. The eigenvalues A±^ and A^ describe the growth 
rates of the perturbations of the velocity parameter a and soliton center x. 

The function (3 1— > sinh(/3) + xdpQpi®) is positive valued. Therefore the real part of the 
eigenvalue A^ is positive and the stationary point labelled — is unstable. Besides, the real 
parts of the eigenvalues A^ and A^ are non-positive for any T > 0. The real parts of A® 
and A^ are also non-positive, which implies that the stationary point labelled + is stable. 

When x is large (say equal to 1), it is important to choose a suitable value to so that 
the soliton parameter (3 + denned by (1421) is small (more exactly, smaller than 0.5). Indeed, 
the theoretical analysis based on the perturbed 1ST is valid only in this case. Furthermore, 
numerical simulations show that (|42p is not a stationary point for larger values of /?+, which 
shows the fundamental limitation in the perturbed 1ST. Within this region of parameters, 
the soliton parameters (3 and a oscillate with the frequency Q + given by f|47|) and they also 
experience an exponential decay with the rate T. The oscillation period and damping rate 
are confirmed by numerical simulations (Fig. [TT]) . 

Moreover, as in the AL case, the propagation of moving solitons is not supported, as 
the soliton velocity decays exponentially to 0. If we denote by «o the initial value of the 
parameter a, then the input soliton converges to its stationary form centered at 




(49) 
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n 

FIG. 12: Here \ = 1, u = 0.3, 5 = 0.022, 7 = 0.02. The initial condition is a soliton with a = cr+, 
x+ = 0, a = 0.1, and j3 = /3+ = 0.45. We plot the soliton profiles |<in(i)| at different times, which 
exhibits the trapping of the moving soliton. The solid line is the theoretical stable stationary 
soliton centered at xf = 4.82 given by (|4U|h 

Notice that we have neglected higher-order term (in j3) in this expression, which is consistent 
with the previous hypotheses and which gives a very accurate prediction for the final soliton 
position (see Fig. IT21). 

V. CONCLUSION 

In this paper we have have investigated the existence and stability properties of new types 
of bright discrete solitons in discrete nonlinear Schrodinger-type models with damping and 
strong rapid drive. Stable stationary solitons are exhibited in the case of a general cubic 
nonlinearity. If the nonlinearity has the special AL form, then stationary solitons, moving 
solitons, and periodic trains of solitons are found to be stable solutions of the system. 
These results have been obtained by applying a perturbed inverse scattering transform 
to the averaged equation and confirmed by numerical simulations. This means that the 
inverse scattering theory is useful for probing the parameter space and exhibiting interesting 
phenomena. One of the problems that should be addressed for future consideration is the 
mechanisms responsible for the instablities of the periodic cn and dn solutions for large drive, 
which seem to be different since a chaotic instability appears first in the dn case, while new 
patterns with different periods appear in the cn case. 
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